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Abstract. The epitaxial growth process of a high symmetry surface occurs 
because adatoms meet and nucleate new islands, that eventually coalesce 
and complete atomic layers. During multilayer growth, nucleation usually 
takes place on top of terraces where the geometry of the diffusion process 
is well defined: We have studied in detail the spatiotemporal distribution of 
£T) | nucleation events and the resulting nucleation rate, a quantity of primary 

importance to model experimental results and evaluate diffusion barriers at 
step-edges. We provide rigorous results for irreversible nucleation and we 
assess the limits of mean-field theory (MFT): we show that MFT overes- 
timates the correct result by a factor proportional to the number of times 
an adatom diffusing on the terrace visits an already visited lattice site. In 
this report we aim at giving a simple physical account of our results. 
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1. Introduction 

Crystal growth by Molecular Beam Epitaxy can be schematically depicted [1, 

2] as a process of uniform deposition of particles, their surface diffusion, and 

their incorporation at the growing surface. Incorporation means attachment 

to another diffusing adatom or to a growing island or to a preexisting step: 

The two first cases are typical of a high symmetry surface where steps are 
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continuously created and destroyed by the growth process, while attach- 
ment to a preexisting train of steps is the growth mode of a vicinal surface. 

Here we investigate the process of attachment between diffusing adatoms, 
called 'nucleation' because the so formed dimer may be the nucleus of a new 
terrace [3, 4]. On a high symmetry surface nucleation may be considered 
as the first step of the growth process, followed by the capture of other 
adatoms and by the coalescence of neighbouring islands. This is surely true 
for the very first stages of growth (the so called submonolayer regime) 
where it is possible to separate the three mentioned steps by monitoring 
the adatom and island densities [5]. In this regime nucleation takes place 
on the flat substrate and nucleation events are not independent processes: 
The formation of a nucleus reduces the adatom density in its surroundings 
and therefore disfavors further nucleations. In a sense, nucleation sites repel 
each other. 

Beyond the submonolayer regime, nucleation mainly occurs on top ter- 
races: in this case, therefore, the diffusion process and the possible nucleus 
formation take place in a confined region of well defined geometry. This 
makes the theoretical study of the nucleation process on top of a terrace 
easier than on a flat surface. However, the spatio-temporal distribution of 
nucleation events cannot directly be related to experiments; in order to 
obtain experimentally relevant information it is necessary to complement 
the results for a top terrace with the growth dynamics of such a terrace. 
Because of this the full study of 'terrace nucleation' is not a 'single terrace' 
problem (we will come back to this issue in the next Section). 

We have not mentioned so far the possibility that adatom attachment 
may not be the end of the story. In fact, depending on the substrate temper- 
ature, the intensity of the incoming flux and the surface symmetry, dimer 
formation may be an irreversible process or not. Irreversible nucleation 
means that once two adatoms meet they stop diffusing and do not detach: 
in the following we are going to study this case. Furthermore islands will 
be supposed to be compact: this requires that attachment of adatoms to 
islands is followed by step-edge diffusion allowing the search of high coor- 
dination sites along the edge of the terrace [6] . 

A short report on the results that will be presented here has already ap- 
peared [7] and a paper with all the mathematical details of the calculations 
is being prepared [8]. Here we aim at giving a simple physical description 
of our approach, of the results and of the open questions. We will consider 
mainly two quantities: the nucleation rate to and the spatial distribution 
P(n) of nucleation events. u> is defined as the number of nucleation events 
per unit time on the whole terrace (it is therefore proportional to the area 
L d of the island); P(n) is a normalized quantity (J2 n P{ n ) = 1) that tells 
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how nucleations are spatially distributed on the sites n of the terrace 1 . A 
third quantity is of interest as well [7]: the probability Q(t) that nucleation 
takes place at time t after the deposition of the second adatom. It has 
less experimental relevance than uj and P(n) and therefore it will not be 
discussed here. 

2. Time scales and basic assumptions 

The process of nucleation on a top terrace involves three typical time scales, 
whose detailed discussion can be found in Ref. [9]. If F is the intensity of 
the incoming flux, D the adatom surface diffusion constant, L the linear 
size of the terrace and d = 1,2 the dimension of the terrace, we define 
T~dep = (FL )~ l as the average time interval between deposition events 
on the terrace, Tt r ~ L 2 /D as the typical time taken by the adatom to 
travel through the terrace, and T res as the typical time an adatom stays on 
the terrace before getting off. This last quantity depends on the bound- 
ary conditions at the terrace edge, i.e. on the possible existence of an 
Ehrlich-Schwoebel (ES) barrier pushing back adatoms [10]; by introduc- 
ing the interlayer transport rate D' we can define the so called ES length 
£ ES = (^ — 1) [2] and in terms of it the regimes of weak and strong bar- 
riers correspond to ^ ES <C L and £ ES 3> L, respectively. The time T res 
has to do with the average stationary adatom density p on the island [9]: 
Tres = p/F ~ L(L + ad( E s)/F) where ad is a numerical factor dipending 
on the dimension and the shape of the terrace. It is straightforward that 
Tres ~ Ttr for absorbing boundaries (^ ES = 0). 

The relation Tt r <C Td ep is always verified in realistic growth conditions 
and therefore, according to the value of the ratio £ ES /L we can distinguish 
three different regimes: (i) Ttr ~ T res <C Td ep (weak ES effect); (ii) Ttr "C 
T res *C Tdep (strong ES effect); (Hi) T tr <C Td ep <C T res (infinite ES effect). 

The word 'infinite' for the third regime is to be intended in physical 
terms: it means that if an atom is on the terrace a second atom will surely 
arrive before the first one leaves the terrace, so that the nucleation rate is 
just the inverse of Td ep - 

Since we study irreversible nucleation and Tt r <C Td ep , we can limit 
ourselves to study nucleation as a "two adatoms" process disregarding pro- 
cesses involving three adatoms or more. We can show the validity of this 
assumption in the most unfavorable case of infinite ES barriers: a first atom 
is deposited at time t = and stays on the terrace; after an average time 
Tdep a second atom comes: they meet in a typical time Ttr an d the proba- 
bility that a third atom arrives in the meanwhile is negligible just because 

Ttr I Tdep < 1- 

1 In d — 2 the index n should be meant as a pair of integer indices (n x , n y ). 
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A further issue deserves to be discussed before the illustration of our 
approach. We study quantities concerning a terrace of fixed size L: it is 
obvious that a given island grows in time, but the nucleation rate u(L) is 
evaluated keeping L constant. The expressions for to and P(n) are general 
and do not depend on the details of the growth process that determine the 
actual time dependence of the terrace size L(t). These details do enter in the 
problem if — for example — we want to compute the probability f(t) that 
a nucleation event has occurred before time t [12], because the evaluation 
of f(t) = 1 — exp{— J (Itu)(L(t))} requires the knowledge of L(t). 

3. The nucleation rate: method and results 

We have argued that it is sufficient to study two adatoms processes. This 
means that the following picture applies: Adatoms arrive on the terrace 
at a rate FL d = tT and stay there an average time r res ; a nucleation 
event takes place if an adatom is still on the terrace when the next one is 
coming and they meet before getting out. It is therefore possible to define 
a nucleation probability per atom p nuc such that u> = tT ■ p nuc . If P dcp (t) = 

r <fep ex P( — t/ T dep) is the probability that a second atom is deposited a time 
t after a first one, we have 

roo 

Pnuc = / dtP dcp (t)p nuc (t) , (1) 

Jo 

where p nuc (t) is the probability that the first atom A deposited at time zero 
and the second atom B deposited a time t later meet. 

In the same manner, if we are looking for the spatial probability distri- 
bution P(n) of nucleation sites we need to evaluate the integral 

/•oo 

P(n)= / dtP dep (t)P{n;t) , (2) 

Jo 

where P(n; t) is the distribution evaluated for the first atom A deposited 
at time zero and the second atom B deposited a time t later. 

Both p nuc (t) and P(n; t) depend linearly on the initial 'probability distri- 
butions' for atoms A and B. Therefore we only need to evaluate p nuc (t) and 
P(n;t) for two particles deposited simultaneously with an effective initial 
distribution [7] for particle A: 



pf(n)= / dtP dep (t)p A {n,t) , (3) 

Jo 

where P\{n, t) is the distribution of a particle evolving alone on the terrace 
for time t: it is the dynamical evolution of the single particle probability 
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p A (n,0) = 1/L via the discrete diffusion equation: 

p A (n,£ + l) = ^5> A (n + <5,£) , (4) 

5 

where "n + 5" labels the neighboring sites of "n" . 

We can explain the meaning of Eq. (3) easily: atom A is deposited 
uniformly at time zero; atom B has a probability P dcp (t) to be deposited 
at time t and therefore it has probability P dcp (t) to find atom A on the 
terrace with distribution p A (n,t). It is worth stressing that p A (n,t) is not 
normalized: its sum on all sites, S(t), is the probability that atom A is 
still on the terrace when B comes in. In conclusion, we have to study the 
problem of two atoms diffusing on the terrace, whose starting distributions 
are p a A {n) and p B (n) = 1/L d , respectively. 

The effective distribution carries two pieces of information: its integral 
(I = J2 n PA B ( n )) determines the scaling of the nucleation rate to while its 
normalized version (p s (n)) influences the shape of P(n). 

Using the definition (3) and the explicit knowledge of -P dcp (£) we can 
write down a differential equation for p B A (n), whose solution is possible [8] 
for any i ES . We propose here a simpler argument: The normalization factor 
of p c A {n) is given by X = r~[ e f °°d£exp(— t/r c i ep )S{t) where S(t) is exactly 
the probability that atom A is still on the terrace at time t. For weak ES 
barriers the exponential can be taken as a constant on the time scale of the 
decay of S(t) and 2 = t"7 J °° dtS(t) = T res /Tdep- For strong ES barriers 
S(t) decays exponentially [13] and I = r res /(r res + r^ ep ). The solution of 
the differential equation mentioned above confirms [8] this result for any 

It is therefore possible to write p c A (n) = Tr ff r — p s (n): The question is 
now how p s (n) looks like. In the limit of weak ES barriers we can repeat the 
above argument and find that p B A (n) = T7$^°p A {n,t). The time integral 
of the single particle probability distribution is nothing but the solution 
of the stationary diffusion equation in the presence of a constant flux [7], 
which is known to have a parabolic form [14] . We can conclude that for any 

^ES : 

p f (n) = T ™ Ps (n) , (5) 

* res ~r T dep 

where p s (n) is the normalized steady state distribution. Consequently the 
nucleation probability per atom takes the form: 

Pnuc = J 5 W (6) 

T~res ~r T^ep 

where W is the probability that two adatoms A and B, deposited simul- 
taneously with distributions p A (n) = p s {n) and p B {n) = l/L d , meet before 
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descending. W is almost independent from the exact spatial profiles of 
atoms A and B and it mainly depends on I-es/L and on the space dimen- 
sionality d. 

Let us now assume T res <C Tdep, i- e. consider weak and strong ES barriers 
(for the 'infinite ES barrier' regime see below). Then, 

2W =* —W ~ pL d W ~ piV dis , (7) 

where we have used the relation p = Fr res and the relation [8] W — N dis /L 
between the probability W that two adatoms meet and the number N dis of 
distinct sites visited by an adatom during its diffusion on the terrace. 

The relation p nuc ~ pN dis can be intuitively justified [15] with the follow- 
ing argument: the nucleation probability per atom is given by the number 
of distinct sites (N dis ) visited by each atom times the probability that a 
given site is occupied (p) . This argument breaks down if the average num- 
ber of adatoms present at the same time on the terrace (and equal to pL d ) 
is larger than one; such condition is equivalent to T res 3> Td ep (regime (Hi)). 

The nucleation rate is finally written as: 

W (L) ~ T ^L ~ FL d pN dis (8) 

dep 

and we want to compare it with the well known mean field result [3, 12]: 

w„f = DL d p 2 = FL d pN M , (9) 

where we have made use of the relation [15] p = -pN M , iV all being the total 
number of sites visited by an adatom during its diffusion on the terrace. The 
comparison of Eqs. (8) and (9) is fully transparent: uj mf /uj ~ N M /N dis = A/", 
that is mean field theory overestimates the correct nucleation rate by a 
factor proportional to the number of times an adatom diffusing on the 
island visits an already visited site. 

Our comparison gets complete once we introduce a model where dif- 
fusing adatoms do not interact: even if they meet each adatom keeps dif- 
fusing until they get off. If the average number of fictitious nucleations 
between non interacting adatoms is W m we simply have iV all ~ W m L d and 
w MF = (T res W N i/T$ ), i.e. mean field theory — as expected by the relation 
Umf = DL d p 2 — treats adatoms as independently diffusing particles. 

The above results are valid for r res <?C Td ep because they derive from 
Eq. (7) rather than from the more general Eq. (6). What does it happen 
in the regime of 'infinite' ES effect? In this case W = 1 and p nuc = 1 so 
that u = t^ = FL d = FN dis while u MF = DL d p 2 ~ FL d pN M . Hence 
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uj mf /lo ~ pL d N M /N diB . We can sum up our results for the nucleation rate 
by reporting the ratio oj mf /uj for the three regimes introduced in Sec. 2, 
both in d = 1 and d = 2: this is done in Table 1 where we use the already 
mentioned relation between 7V all and p, the well known results [16] that in 
absence of ES barriers N diB ~ L in d = 1 and iV dis ~ L 2 / In L in d = 2 and 
finally that N dis = L d for £ ES > L. 



TABLE 1. Ratio Wmf/w for the three relevant regimes. For 

(i) and (it), umf/u = N = N^/N^; for (m) wmf/u — pL d N 
where p ~ -^L£es is the average adatom density. 

(i) (ii) (m) 

d=l L fe s F(L£ ES f/D 

d = 2 InL £ E s/i F(L£ ES ) 2 /D 



The result lo mf /u = i ES /L — valid in the two dimensional strong bar- 
rier regime — had already been found in Ref. [9]. In that paper logarith- 
mic corrections were neglected, which corresponds to disregarding the In- 
dependence of the factor W in the weak barriers regime. 

4. The spatial distribution of nucleation events 

In this Section we study what sites are the most favored for nucleation. 
The result is strongly dependent on the exact spatial profiles of the initial 
probability distributions for particles A and B (p A]B (n)): we have justified 
below Eq. (3) that A is distributed as the stationary solution of the diffusion 
equation (p\(n) = Ps{n)) while B is uniformly distributed (p B (n) = 1/L d ). 

The dynamical evolution of two diffusing particles on a d— dimensional 
terrace can be easily mapped on the problem of a single walker in a space 
of dimension d' = 2d. If m and n label the positions of the two atoms a 
nucleation event occurs when m = n, i.e. when the "single walker" crosses 
the diagonal of a square terrace if d = 1 {d! = 2) or the 'diagonal plane' of 
a four dimensional hypercube if d = 2 {d! = 4). 

In d = 1 we have solved analytically [7] the problem for the limit cases 
of zero and infinite ES effect; for any value of £ ES we can compute exact 
numerical results for P(n) both in d = 1 and d = 2 [8]. Our main results 
in d = 2 are reproduced in Fig. 1 where we give the resulting spatial 
distribution of nucleation events along the diagonal of a square terrace 
and we compare it with Mean Field Theory: P(n) ~ Ps(n). In absence of 
ES barriers (Fig. la) MFT works remarkably well but its accuracy gets 
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Figure 1. P(ri) along the diagonal of a square terrace of size L = 20 for (a) £es = 0, 
(b) £es = 4, (c) £es = 20, (d) ^es = 100. Full diamonds: exact theory; open squares: 
mean field theory; full line: "independent adatoms" model. 



worse with increasing £ ES (Figs. lb,c) and MFT fails completely for large 
ES barriers (Fig. Id): in this case MFT predicts that P(n) gets flat while 
a clearly rounded shape is obtained. 

Fig. 1 confirms that MFT is equivalent to an "independent adatoms" 
model and this is the reason of its failure with increasing £ ES : f° r 4s ^ L 
independent adatoms perform many more fictitious nucleations than for 
4s < L (see Table 1). 



5. Conclusions 

We think that two main aspects of our work should be emphasized. Limits 
and inaccuracies of mean field theory have been pointed out by several au- 
thors during the years [17, 11, 9, 18, 19] but we have performed a rigorous 
study for any value of the ES barriers in one and two dimensions. Further- 
more we study the spatial distribution P(n) of nucleation events and we 
are not aware of previous similar analyses. The second aspect is that we 
provide a simple physical interpretation of the failure of mean field theory: 
it counts all the nucleations of two independently diffusing adatoms. 

The nucleation rate is of primary experimental relevance because most 
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of the methods to determine step-edge barriers — or to evaluate other quan- 
tities derived from £ ES — require the knowledge of u>(L) [12, 20]. It would 
be therefore useful to reconsider such derivations which made use of uj mf : 
this has already been started in Ref. [9]. 

As for P(n), a direct experimental determination is extremely com- 
plicated [21]. Nonetheless its knowledge has a theoretical interest because 
P(n) may enter in mesoscopic models of crystal growth [22, 23], i.e. in 
models where surface diffusion is not taken into account microscopically 
but through a mesoscopic surface current plus the rule for nucleating new 
terraces. 

Finally we mention a few extensions of our work that are presently 
in progress: Firstly, the problem of nucleation between adatoms of differ- 
ent species, i.e. particles having different diffusion constants. Secondly, the 
study of P(n) in the limit L — > oo in order to understand what features are 
maintained in this limit. Finally, the nucleation on top of a fractal terrace. 
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